Power Dependent Lineshape Corrections for Quantitative Spectroscopy 
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The Voigt profile - a convolution of a Gaussian and a Lorentzian - accurately describes the ab- 
sorption fines of atomic and molecular gases at low probe powers. Fitting such to experimental 
spectra yields botfi tfie Lorentzian natural linewidth and the Gaussian Doppler broadening. How- 
ever, as the probe power increases saturation effects introduce spurious power dependence into the 
fitted Doppler width. Using a simple atomic model, we calculate power-dependent corrections to 
the Voigt profile, which are parametrized by the Gaussian Doppler width, the Lorentzian natu- 
ral linewidth, and the optical depth. We show numerically and experimentally that including the 
correction term substantially reduces the spurious power dependence in the fitted Gaussian width. 



Vapour cell spectroscopy is important for determin- 
ing the properties of atomic or molecular transitions 
[Tl[2], which are well described by the Voigt profile, the 
workhorse of transmission spectroscopy. Experimental 
progress has reached the point that intensity dependent 
corrections to the Voigt profile are becoming significant. 
For example, molecular fingerprinting and other trace gas 
measurement applications [3l |4] using direct frequency 
comb spectroscopy (DFCS) 3] require a quantitative link 
between absorption depths and gas concentrations, par- 
ticularly if the probe laser intensity levels are increased 
to improve signal-to-noise in challenging environments. 
While absorption sensitivities measured using DFCS can 
be determined with quantum- noise limited precision [6], 
the effects of optical pumping on the observed depths 
and line-shapes remain to be quantified. Similarly, op- 
tical pumping can limit accurate determinations of the 
Boltzmann constant, fcs, from measurements of Doppler- 
broadened line shapes [71IT3]. To optimise the trade-off 
between signal-to-noise limitations at very low powers, 
and systematic optical pumping effects at higher powers, 
an accurate theory of the line shape dependence on probe 
intensity is needed. 

In the limit of low probe powers, where atomic popu- 
lations are hardly perturbed from their thermodynamic 
equilibrium values, transmission spectra are well de- 
scribed by an exponentiated Voigt function [TH [T5] . 
7iin(^) = g-azV,,{A)^ where a is the absorptivity, z is 
the length of the vapour cell, A is the detuning from res- 
onance (note: all times and frequencies are expressed in 
units of F = Fnaturai/2), v = "f/T is the non-dimensional 
Doppler width, and y is a convolution of a Gaussian and 
a Lorentzian: 
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transmission spectrum very accurately; the Gaussian 
component arises from the Doppler shifts due to the 
atomic Maxwell-Boltzmann velocity distribution, while 
the Lorentzian relates to atomic relaxation processes. 
Fitting Tiin to experimental measurements yields 7, F 
and the optical depth. However, as the probe power 
increases, perturbations to the equilibrium atomic pop- 
ulation distribution become significant and Tiin fails to 
accurately represent the transmission spectrum. In the 
limit of large probe powers, atomic populations depart 
from thermodynamic equilibrium: in two level atoms the 
ground and excited states tend to equalise; in three level 
systems population may be transferred to an optically 
inactive state. In either case, fitting Tiin to the measured 
spectrum yields incorrect values for 7 and F, which ac- 
quire a spurious intensity dependence. 

In this letter, we derive the intensity dependent cor- 
rections to Tiin, which can be computed perturbatively 
in powers of the laser intensity. We then show, both 
theoretically and experimentally, that fitting to the cor- 
rected form yields Doppler widths that are independent 
of intensity as they should be. This result will enable ac- 
curate measurements of Doppler broadening in quantita- 
tive spectroscopy at much higher probe intensities, where 
saturation effects are not negligible. This in turn greatly 
enhances signal-to-noise in a given integration time. 

We begin by computing the spectral dependence on 
the atomic populations for a three level atom, consisting 
of two ground states: one optically active, |1), and the 
other optically inactive, |2), and an excited state |3), as 
discussed in ^16,. Transitions between states |1) and |3) 
are optically driven, and state |3) can relax to either of 
the ground states. The population rate equations are 
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(1) 



where erfc(z) = 2 J e * dtj ^. 

At low probe powers, Tiin describes the measured 



where P = {Pi,P2,P3}, fl is the atomic Rabi frequency 
(proportional to the electric field amplitude), (3 is the 




FIG. 1. (Colour online) Typical forni for the Voigt function 
(dark, blue), and the correction term multiplying q in Eq. (Im 
(light, red), for v = 10, p = 20. 
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FIG. 2. (Colour online) Fits to numerically simulated spectra, 
showing the ratio of fitted Doppler to lorentz widths demon- 
strate the improved performance of fits using the correction 
term. 



branching ratio from state |3) to 1 1) and A is the detuning 
between the laser frequency and the atomic transition. 

The light is absorbed as it propagates through the 
atomic medium, and the axial evolution of the field is 
governed by ^16, 



dlnVL^ 



dz 



= 2a 



dA. / P(A„J. (2) 



Here P(A„J = P3 ( A^ J - Pi ( A„ J and the notation A„^ 
indicates that a given atoms sees a detuning which de- 
pends on its axial velocity via the Doppler shift. If the 
field intensity is weak, so that on resonance populations 
are negligibly perturbed from thermal equilibrium, then 
P = -1/2, and the integral yields n^{z) = flge"^^-^'^). 
Since T(z) = i}'^{z)/Ql, we recover T = Tiin- 

As n^ grows, perturbations to P become significant, 
and we should compute corrections to P using Eq. (fl]) . In 
the simplest case where il^ is time-independent Eq. (Il]) 
can be solved analytically. The expression is cumber- 
some, however it may be straightforwardly expanded in 
powers of fl^ to yield 



P{A,J = -l+q{t) 
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where q{t) w (1 + /3 + 2 t(l - /3))/8, and t is the time for 
which the atom is exposed to the field. Substituting this 
result into Eq. (pi) yields 



dz 
where 



= -aK(A) + aq{t) ^^^^(A) Q^ + OiQ^) (4) 



FJ")(A) = 
is a generalisation of the Voigt profile 
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V'^-^^ may be evaluated by noting that 
= lim 
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(l + a;2)2 a^i [ (a2 - 1) (a2a;2 + 1) (1 - a2) (a;2 + 1) 

Convolving this sum of two Lorcntzians with a Gaussian 
thus yields a sum of Voigt functions. 
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Higher order correction terms can be calculated itera- 
tively, using the same method. 

Equation (B can be solved analytically to give 



T = 



zV„{A) 
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Expanding this expression, and recalling that all quan- 
tities are nondimensionalised in units of T we explicitly 
include the Lorentz width as a parameter 

J- _ gpV„(A/r)+q(e'=^-<^/r')_i) yJ2)(A/r)/K(A/r) ^-g^ 

where p — —az and q = —qflQ. This correction, lin- 
ear in g oc /, is the central result of this Letter. If ei- 
ther the intensity (ex q) is small, or if the optical depth 
(ex p) is small, the correction vanishes and the well-known 
Voigt profile is recovered. In the latter case, saturation 
effects may be substantial, however the optical depth of 
the sample is sufficiently short that only a small frac- 
tion of incident photons are absorbed. In what follows 
we demonstrate both numerically and experimentally the 
advantages of using Eq. ([6]) when fitting spectra for which 
saturation effects become significant. 
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FIG. 3. (Colour online) Fits to (left) Simulated data using Eq. pi, with /? = 1/2, f = 20, ;/ = 100; (right) experimental data 
taken from [^. (top) Fitted Doppler widths, extracted from fits to Voigt only or Voigt plus correction, with fixed or free 
Lorentzian widths, (bottom) Fitted Lorentzian widths, normalised to the natural linewidth. 



We note that the derivation of the correction term as- 
sumed the Rabi frequency was time independent. How- 
ever, in an experiment atoms crossing through a probe 
laser beam see a time-dependent field. Nevertheless, 
Eq. (Is]) and (|4| are still correct, albeit with some more 
complicated time dependence in q{t\ and taking f2 to 
represent a typical scale of the local Rabi frequency (e.g. 
that corresponding to the peak intensity in the probe). 
It follows that the form of the correction term in Eq. (l6]) 
is also valid, where the complicated time dependence in 
q is simply absorbed into the parameter q. 

A typical Voigt profile and the correction are shown in 
Fig. [l] In fitting to spectra, v, p and q are treated as 
fitting parameters. We also numerically investigate leav- 
ing F free, which heuristically accounts for higher order 
intensity dependences that we have otherwise neglected; 
variation in F implies a breakdown in the fitting model. 

To evaluate the efficacy of the correction, we numer- 
ically simulate an absorption spectrum using Eq. (pi) 
(keeping terms in P up to 0{V,^)), and then attempt 
to extract the simulation parameters by fitting Eq. (pi). 
For comparison, we compare fits with Voigt only (q = 0) 
and Voigt plus correction (g free); in combination with 
a fixed Lorentzian (F = 0) and free Lorentzian (F free). 



Fig. p] shows the deviation of Vfn = 7fit/r from the cor- 
rect value, I/, as a function of intensity (ex (ilo/r)^), for 
the different fitting forms. As expected, all fitting forms 
yield the correct value for i/fn as ilg ~^ 0- Without the 
correction (■ and •), the fit converges to the correct 
value as Hq. This improves to Q,q when the correction is 
included (♦ and A). 

The fitting performance is further illustrated in 
Fig. p] (left), which shows the dependence of the fitted 
parameters Vfn/i^ (top) and F (bottom) as a function 
of the probe intensity (ex (JIq/T)'^, for a particular rep- 
resentative choice of parameters (/3 = 1/2, t = 20, 
h' — 100). Also shown is the fitted Lorentzian component. 
Fig. p^ (bottom), demonstrating apparent broadening of 
the atomic lifetime, even for intensities substantially be- 
low the saturation intensity. This phenomenon has been 
observed experimentally [T7] . 

Qualitatively, at low powers the uncorrected, Voigt- 
only fits (I and •) systematically underestimate v over 
a range of laser powers, and demonstrate a pronounced 
minimum, as illustrated in Fig. pTtop) . In the case where 
the Lorentzian is held fixed (•), the Voigt-only fit may 
over- or under- estimate the Doppler width, depending 
on the parameters, i.e. the initial slope may be positive 



or negative. For the parameter choices in the simulated 
spectra, Fig. [s] (left), the initial slope is slightly negative, 
and there is a noticeable minimum in Vf^t/v- 

We now demonstrate the effectiveness of using the 
corrected transmission function for extracting the true 
atomic velocity distribution from experimental data. 
Spectra were recorded in a 7 cm long Caesium vapour 
cell at approximately 295 K at different probe powers, 
using a coUimated 2 mm diameter laser. At each power 
setting, 16 spectral scans were taken, and each scan was 
fitted using Eq. ^ to yields estimates for 7fit and Ffit. 
For each scan, the temperature was recorded using a plat- 
inum resistance thermometer (PRT), accurate to ±100 
ppm. Each fitted Doppler width, 7fit, was normalised 
by the width expected given the temperature measured 
using the PRT, 7prt = Wprobe\/2A;BTpRT/TO/c. 

To qualitatively compare the fitting results from the 
simulated spectra with the experimental results, we con- 
vert the laser power to an average intensity by dividing 
the power by the beam area, 7r(l mm)^, and then nor- 
malise this average intensity by the saturation intensity 
^sat = 2.51 mW/cm [TS], so that the nondimensional ex- 
perimental parameter, 2///sat, is equivalent to the nondi- 
mensional simulation parameter (O/F)^. 

The results are shown in Fig. [s] (right). For 
2^/^sat < 0.1, the extracted Doppler widths are substan- 
tially better when using the correction than not, as shown 
in Fig. p^ (right, top). This is in qualitative agreement 
with the simulated data. For the lowest powers used, the 
corrected estimates for the Doppler width are at least 
an order of magnitude better than the uncorrected form, 
limited by the accuracy of the PRT calibration, and sta- 
tistical uncertainty from the small sample size. 

Comparing Fig. [s] (left) and (right) shows that the sim- 
ulated results are in qualitative agreement with the ex- 
perimentally derived results in other respects as well, in- 
cluding the relative locations of minima, slopes and cross- 
ing points. Quantitative differences arise from the simpli- 
fied physical model underpinning the simulated spectra, 
in which the distribution of atomic crossing times and 
intensity profiles is replaced by the average crossing time 
and the average intensity. While not unphysical, this 
differs from the experimental conditions. 

In thermometry, where the objective is to extract very 
accurate Doppler widths from spectral data, we see from 
Fig. [2] that it is important to use the Voigt correction. 
For instance, if 2I/Isa,t ~ 0.003 then a fit to a Voigt- 
only profile has a systematic bias in the fitted Doppler 
width at the level of several hundred ppm. In contrast, 
using the corrected form yields an estimate of the Doppler 
width that is accurate to ^ 1 ppm, which is comparable 
to the current thermometric state-of-the-art. 

It follows that under the comparable experimental con- 
ditions represented in Fig. IS] (right), to reduce the sys- 
tematic bias in the fitted Doppler width to 1 ppm (when 
utilising the corrected form) , the probe beam should op- 



erate at 2///sat — 0.003, corresponding to a laser power 
of 100 nW in a 2 mm diameter beam. 

One important application of Doppler spectroscopy 
is to contribute to the CODATA redefinition of Boltz- 
mann's constant [H [SI HH US] , in which the vapour cell 
is at equilibrium with a triple-point-of-water reference at 
273.16 K. At this temperature the optical depth of the 
cell is ~ 10 times smaller than under the experimental 
conditions reported here. As such, the intensity may be 
increased by the same factor, to ~ 1 /iW, while still main- 
taining the systematic Doppler error at ~ 1 ppm. 

In conclusion, we have derived the power-dependent 
correction to the Voigt profile, and demonstrated numer- 
ically and experimentally that including the correction 
term yields much better estimates of the width of the 
underlying Gaussian process, substantially reducing the 
spurious power dependences that arise from an uncor- 
rected Voigt profile. We anticipate that this will find di- 
rect application in the near-term to high-precision spec- 
troscopy and thermometry of atomic vapour cells. 
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